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Abstract 

An eigenvalue problem relevant for non-linear sigma model with singular metric is considered. We 
prove the existence of a non-degenerate pure point spectrum for all finite values of the size R of the 
system. In the infrared (IR) regime (large R) the eigenvalues admit a power series expansion around 
IR critical point R — > oo. We compute high order coefficients and prove that the series converges 
for all finite values of R. In the ultraviolet (UV) limit the spectrum condenses into a continuum 
spectrum with a set of residual bound states. The spectrum agrees nicely with the central charge 
computed by the Thermodynamic Bethc Ansatz method. 
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1. Introduction 

The non-linear sigma models in two-dimensional (2D) space time are widely used in field 
theory as continuous models of two-dimensional spin systems (see e.g. Ref.s. ^ 121111111]) as 
well as in relation to string theory (e.g. Ref.s. [HI El 13 El ) ■ The general 2D sigma model 
(SM) is defined through the action 

(1) A[G] = l -f GijiX) d^X 1 d^Xi d 2 x 

where the coordinates x^, fj, = 1, 2 span a 2D flat space-time, while the fields X 1 , i = 1, d 
are coordinates in a (i-dimensional Riemannian manifold called target space. The symmetric 
matrix Gij is the corresponding metric tensor. 

The standard approach to 2D sigma models is perturbation theory. If the curvature of 
Gij is small, one can use the following one-loop renormalization group equation from Ref. [3] 



(2! 



where t is the RG "time" (the logarithm of scale) and Rij is the Ricci tensor of G. 

The simplest examples of 2D sigma models are SM with two-dimensional target space 
(d = 2). In this case Rij = H8ij where 71 is the scalar curvature. Then we can always 
choose (at least locally) conformal coordinates for which 

(3) G^ = e* 5ij 

with a single function <£. 

An important role in the analysis of 2D sigma models is played by the effective central 
charge c{R). This dimensionless function contains the information about the UV and IR 
properties of the theory and it is related to the ground state energy Eq(R) of the corre- 
sponding quantum system, living on a finite space circle of length R: 



(4) *W>~^. 

For SM with two-dimensional target space cjjv = c (0) = 2. 

In integrable theories this quantity can be calculated exactly using Thermodynamic Bethe 
Ansatz (tba) equations Ref . |§1 ITU] . This problem is however much more complicated for the 
excited levels Ei(R), so it is useful to have some independent approach for their calculation. 
It was shown in Ref. jllj that for the sigma models with d = 2, in one loop approximation 
(EqEJ), this problem can be reduced to the eigenvalue problem for the operator: 

(5) h = -\vl + \n t , /^ = ^% 4 . 

Here Vf = e~^(d 2 + dy) is the Laplace operator and Tit is the scalar curvature in the SM 
metric renormalized at the scale R: 

(6) t-t = logRA 

where Ao is the normalization parameter. This operator is self-conjugate with respect to 
the scalar product in the SM metric: 

(7) (*i,* 2 ) = / *i* a e*da;dy. 
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The effective central charge c(R) in one loop approximation can be expressed through the 
lowest eigenvalue : 

(8) c(R) = 2 - e (R) 

and the excited levels Ei(R) = E$(R) + ir(ei(R) — eo(R))/6R. We note that if in the IR limit 
SM flows to the critical point described by conformal field theory (CFT) then the numbers 
= (ej(oo) — eo(oo))/24 coincide with conformal dimensions of the fields in this CFT. 

The eigenvalue problem (Eq. ©) with the natural scalar product provided by the metric 
Gij(X) can be applied to the analysis of 2D sigma models with target space of arbitrary 
dimension d . It follows from Zamolodchikov's c-theorem that the effective central charge 
defined in Ref.|12| through the correlation functions of the energy momentum tensor is 
non-increasing as a function of the scale R. The effective central charge defined by Eq. (jHJ) 
(with d — eo(R) in the r.h.s.) also satisfies this remarkable property, which follows from one 
of the results of Ref. 13j, where it was shown that the lowest eigenvalue of the operator h 
is a non-decreasing function of R. 

A well known solution to Eq. (j2J) (see E,ef.|llj) defines the axially symmetric metric of 
the "sausage" SM (an integrable deformation of the 0(3) non-linear sigma model), which 
is described by the action: 

W •Assg-J a ( t ) + b(t) cosh2Y a X 

where a(t) = z^coth2u, b(t) = z//sinh2u, and u = z^(to — t)/4n. It is easy to see from the 
explicit form of the metric © that operator hjv depends only on the variable u and does 
not depend on parameter v. It means that 

(io) <•»(*) = ^*r». 

where after the substitution ^ = ex.p(imx)^ m (y) (to £ Z) the scaling function k^ S9 (u) is 
the eigenvalue of the Sturm-Liouville problem: 



(11) 



2 2 1 + cosh 2u cosh 2y ^k^ 9 (u) sinhu 
~ y 



(cosh 2u + cosh 2y) 2 cosh 2u + cosh 2y 



<JrO0 = 



where the eigenfunctions have finite norm according to Eq.(|7j). For the ground state 
Kq S9 (u) this problem was studied in Ref. [TT] . 

In this paper we consider the eigenvalue problem for the sigma model which correspond 
to another solution of RG Eq. (|2*|l. This solution can be obtained by analytic continuation 
Y — > Y + iir/4, u — > u + in/ 4 from the solution for the sausage model. The corresponding 
action can be written as: 



a(t) + P(t) sinh2y 

where a(t) = i/tanh2n, (3{t) = i//cosh2u, and u = v(to — t)/4-K. This metric has a 
singularity at Y = —u. It means that the coordinate Y in target space should be considered 
only in the region Y > —u. The curvature 1Z also has singularity at this point. However for 
small values of parameter v the curvature is not small only in the narrow region {SY ~ v) in 
the vicinity of singularity. A more careful analysis shows that the one loop approximation 
is valid for the calculation of the observables in SM (|12j) . The relative correction to one 
loop approximation as well as in the sausage SM has the order i/log(l/V). The eigenvalue 
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equation for the scaling functions Ki(u) (ei(R) 



UKi{u) / '47r) has now the form: 



(13) 




1 — sinh 2u sinh 2y g- K% (u) cosh u 




(sinh 2u + sinh 2y) 2 sinh 2u + sinh 2y 



The solution (y) should satisfy now the boundary condition 



(14) 



(y + u) 



1/2 



fi- 



ll 



and it must be square integrable with respect to the natural norm 



(15) 




'OO 




dy . 



— a. 



sinh 2u + sinh 2y 



In the IR limit u —* — oo the metric ()12j) has an asymptotic which can be written in the 
form of Eq. (J3J) with exp{— &ir} = |z/(exp{2Z} — 1), with Z = Y + u. For discrete values 
of the parameter v = 4tt/N the SM with this metric can be derived from the SU(2) level 
N WZW models by gauging U(l) symmetry (see Ref. |141 115) for details). The resulting 
SU(2)n /U(l) coset model described by the SM with metric exp($>m)5ij coincides with Zjv 
parafermionic CFT of E,ef.|16j. For general values of u the quantum field theory (QFT) 
corresponding to SM 1|12|) can be considered as the deformation of parafermionic CFT. It is 
natural to expect that it will be a massless theory describing the RG flow from rather non- 
trivial UV field theory (which is also well defined for the same discrete values of parameter 
v) with cjjv = 2 to the parafermionic CFT with cjr = 2 — 6/(N + 2) in the IR limit. 
The scaling functions Ki(u) in this case describe the RG dynamics of energy levels from UV 
regime to IR asymptotics, where they define (with relative accuracy 0(1/N)) the spectrum 
of anomalous dimensions of parafermionic CFT (see Appendix C). 

This massless QFT is integrable and can be described in terms of factorized scattering 
theory for massless excitations. For v = At: /N with N > 3 the ground state energy Eq.Q of 
the SM (|12ft (as well as that of the sausage SM Q, of Ref.fllJ admits an exact calculation 
by tba method. These equations will be described in Sec. 03 

Both eigenvalue problems Eq.(ll,13) are believed to have a purely discrete spectrum. 
The ground state eigenvalue Kq S9 of Eq. (JTT|) was studied in Ref. , where the asymptotics 
of this function in the regimes u — * and u — > oo were found. This function was studied 
numerically in Ref. ^7] and the result was in perfect agreement with the scaling function 
calculated from tba equations (see Sec. 5). Actually it is rather easy today to attack the 
spectral problem Eq. (|11|) by using sparse matrix techniques available in mathematical 
libraries. Good accuracy can be achieved by introducing a multi-grid discretization and 
using Richardson extrapolation. This approach, however, is not immediately applicable to 
new equation 1)13(1 due to its singular nature. In the following we shall bring the equation 
to a form which is suitable for a detailed perturbative analysis (Sec. HJ), to a second form 
which allows an accurate asymptotic analysis in the UV regime (Lame form, SecOJ), as well 
to a third form, more suitable for a purely numerical approach (Sec. 0J). Finally (Sec. EJ) 
we shall exhibit the matching of the ground state of Eq. (|13|) with the central charge of 
the modified TBA system whose structure is given in Fig. 8. The interested reader will find 
some further mathematical details in the Appendices. 
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2. The connection with Heun's equation 

To bring Eq. ()13|) to a more manageable form, we begin by re-absorbing the the integration 
measure into the wave function. By defining 

*G0 = y/pty) <l>(v) , 

(16) p{y) = sinh 2y + sinh 2u 
we find 

(17) -A L)M + ( m 2 p (2/) - smh2y) #y) = ± K cosh2n 0(y) 

Putting x = e~ 2 ^ +n ^ and it; = — e _4n , Eq. Q17JI is transformed into the following 

, , „, „ / 1 1 \ , (1 - m 2 )x 2 - \k (1 - w)x + m 2 + w , s 

(is) ^ ( x) + (— + — ) m + i L-^-L^ «,) = o 

of the Fuchsian type. A further substitution <j> —* yjx f(x) reduces Eq. (|18jl to the form 

(19) rw + (— + + -^-) /'(x) + < 1+ ;f"\ / w = o 

\ x x — l x — wj x[x — l){x — w) 

where the so-called accessory parameter q is given by 

(20) q = Jj(«(l-«;) + 6(l + u;)(l + 2m)) . 

This equation was analyzed by Heun in 1888 (see Ref. ^3E3) wno considered a general 
linear differential equation of the second order with four Fuchsian singularities. In Heun's 
notation, the solution is formally given by F(w; q; 1 + m, 1, 1 + m, 1; x), but this is of little 
use in practice. We gain some insight from the fact that the limit w — > — oo is a case of 
confluence of singularities which takes us back to the hyper geometric equation (see Ref. ^H] 
for a general treatment). 

2.1. The eigenvalue problem in algebraic form. It is well known that series solutions 
for Heun's equations can be most conveniently constructed using the basis of Jacobi poly- 
nomials Pn m '°^(l — 2x) (see e.g. |20j . Vol. III). We are now going to show how to solve the 
eigenvalue problem by exploiting this favorable basis: the problem will reduce to finding 
the spectrum of an (infinite-dimensional) tridiagonal matrix for which efficient algorithms 
are well-known to exist |2j . Let us consider Eq. (|19j) : after setting y = 1 — 2x , we can 

easily expand the solution in a series of Jacobi polynomials V™ = P^"" 1 °^ (y) by converting 
the differential equation in the form 

(21) Hf= j(l-2«;-y)A>(A> + m + l) + (l-y 2 )|--(l + m)y}/ = (2q - m - 1) / 

where NV™ = nV™. Now we can use the basic properties of Jacobi's polynomials (see 
e.g. [13, Vol.11, [22]) to reduce the operator H to the form Hf = (1 - 2w) H f + V / whose 
action on the basis vectors is particularly simple: 

H V™ = n(n + m + l)V™ 

= (2 + m)n 2 + (m + l)(m + 2)n + m(m + 1) 

(22) V ' n (m + 2n)(m + 2n + 2) 

n 2 {m + n) 2 (n + l) 2 (m + n + l) 2 

(m + 2n)(m + 2n + l) " _1 (m + 2n + l)(m + 2n + 2) n+1 

We can now conveniently study the spectrum of q using this tridiagonal matrix represen- 
tation, by applying, for instance, the technique of Sturm sequences and bisection |2J. The 
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matrix representation also lends itself to a very simple perturbation series expansion, as we 
discuss in the next section. We shall have to refer to the matrix representation of H in the 

orthonormal basis = (-) n J ^q^TT^™ as V nn r - 

2 (n + l) 2 (n + m+l) 2 

Vn+l,n - 



V, 



n—l,n 



(2n + m + 2)^(2™ + m + 2) 2 - 1 
2n 2 (n + m) 2 



(2n + m)y / (2n + m) 2 - 1 

diagonal terms being unchanged, but the former V, being rational in its indices, is more 
convenient for the calculation of perturbative coefficients; we show in Appendix A, that we 
can use V without modifying the standard algorithm. 

2.2. Perturbation theory. The tridiagonal matrix representation of H, when rewritten 
as H = e~ 1 (Hq + eV), can be used to calculate a perturbative expansion in the parameter 

(24) e = (1 - 2 w)- 1 = (1 + 2c" 4 ")- 1 . 

The convergence of the expansion is governed by the Kato-Rellich theorem^: let there 
exist constants a, b such that ||V0|| < a\\(f>\\ + 6[|Ho^|| (||-|| denotes Li-norm). Then the 
perturbative expansion defines a regular analytic function for \be\ < 1. For a tridiagonal 
matrix it's not so difficult to find norm estimates; in our case it's simple algebra to check that 
the column sums of the matrix elements of V coincide with the diagonal matrix elements of 
Ho, up to an additive constant, therefore we have 

(25) V = S (Ho + (m + 1)1) 

with S a stochastic matrix" and I the identity matrix. Hence it follows 

(26) ||V^|| = \\S(H + (m + 1)I)0|| < \\HoH + (m + 1) 



which implies that the perturbative series will converge for e < 1, that is for all u. Actually 
we can say more: let us denote by R(H, n) = (H + the resolvent operator; by applying 
Schur's Test to the matrix R(Ho,n) V (the symmetric version of V) one concludes that 

(27) \\R{Ho,li)V\\<l 

if n > m + 13/8. Now the resolvent of H satisfies Lippman-Schwinger equation 

(28) R(H, n) = R(H ,n) - eR(H , n) V R(H, ft) 
which for |e| < 1 can be inverted to give 

(29) R(H, n) = (1 + eR(Ho, /i) V)' 1 R(H , n) . 

Since R(Ho,fJ-) is a compact operator, Eq. ()29() is telling us that the resolvent R(H,/jl) is 
itself a compact operator, being the product of a compact operator with a bounded one. 
This implies that the spectrum of H is purely discrete. 

The expansion can now be computed rather easily by the standard recursive algorithm 
(see Appendix A). Details on the series expansion can be found in Appendix B, where we 
prefer to adopt a different parameter which naturally appears in the Lame formulation of 
next section, namely 

(30) A = (l + e- 4u )- 1 = ^, e = -^r. 
y J K J 1 + e' 2-A 



t See for instance 
Hie. E,5 tJ = 1 
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The expansion in powers of A turns out to be simpler and with better convergence properties; 
indeed the substitution e — > A is just a special case of Euler's (E, q)-method J2Jj|. We present 
just a sample of the infinite number of different series expansions, since we believe that 
nobody would like to copy them from paper but would rather prefer getting the code which 
generated the expansion**. The first few terms for the ground state value of k at fixed to 
are the following 

, . i 2m , 4(m + 1) 3 x9 

31 = l + 2m -X- - \ ; A 2 - 

D to + 2 (to + 2J' 3 (m + 3) 

8(to + 1) 3 (2to 2 + 5to + 4) ^ 4 

(to + 2) 5 (to + 3)(to + 4) 1 ' 

while for the excited states, after putting j = |m + n, we have (see also Appendix C) 



[4j(j + l)-r^' 2 

->>»,„ i -./ + u - /// 8j(JTT) 



■ i ii 777, 

(32) ±K m , n = (2j + 1) 2 -to 2 



1 



2 9 (2j + 1) 



( 4(i + 1)2 _ m 2)4 (4j 2_ m2) 4- ^ ^ 



(i + 1)3(2; + 3) i 3 (2j - 1) 

Using floating point arithmetic we may quickly explore very high orders, with due atten- 
tion to truncation errors which accumulate along the iteration. The asymptotic behavior 
of the coefficients shows very clearly a limit c n +i/c n — > 1, confirming that the series con- 
verges in the unit circle, which means in the domain \e /(l + e 4u )\ < 1. In the complex u 
plane this is a domain which includes the whole real axis. However as we venture along the 
positive real axis, the convergence is critically slowed down: to go deep in the U V region we 
may be obliged to sum a really huge number of terms, or try some resummation, e.g. via 
Pade approximants. Since high order coefficients are easily computed, however, we may try 
to extract the asymptotic behavior of k(u) for large positive u by analyzing the asymptotic 
behavior of the coefficients. For example we can verify that the UV asymptotics of the 
ground state eigenvalue kq(u) coincides with that of function Kq S£i (u) (see SecEJ) and has a 
form: 

3tt 2 

(33) K (u) = + 0(l/u 5 ). 

2(u + log4j z 

When expressed in terms of e this formula can be expanded in a power series and the 
coefficients compared to those coming from perturbation theory. Sub-dominant terms tend 
to mask the simple n -1 (logn) -3 behavior one should expect; we find that it is more accurate 
to compare the Taylor coefficients of the function l/(4q — 1), which is actually diverging 
at e — * 1, with those of its leading term 0(log(l — e) 2 ). However, to make the comparison 
even more transparent, we may look for a special function whose behavior at e — * 1 is the 
simplest possible. Let us observe that 2q — 1 turns out to be an odd function of e; the new 
function yje/ (4q — 2 + e) has the simple 0(log(l — e)) leading singular behavior at e — > 1, 
and it is an even function of e. We argue that its leading behavior should then be 

(34) T[ % — !— log ^ 



2vre V 1 - 

In fact we find that the expansion of T matches the perturbative series with a high accuracy 
(see Fig. 1 where the deviation is magnified 10 5 times). The asymptotic behavior will 
be recovered in a very precise way numerically in Sec. |2J hence its extraction from the 



"Matlab and Mathematica codes are available at the web site www.fis.unipr.it/~onofri. 
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Figure 1 . Relative deviation of T 's expansion coefficients from those of its 
leading asymptotics (xW 5 ) (Eg 131)). 



perturbative series appears of purely academic interest. Anyhow, assuming Eq. 1)3*1)) it 
follows from Eq. ((20)) 

(35) «= — «6(4q-l) = 6/T 

since |log((l + e)/(l - e)) = ±log(l + e 4u ) « it as e -> 1. 

2.3. Sausage model equation. The sausage model was our starting point. Now we go 
back to it and show how it fits into the correspondence with Heun's equation. The scaling 
function is defined in Eq. 1)11)1 . We shall now look for its algebraic equivalent as we did in 
Sec. Since the two equations are related by analytic continuation, it will not come out as 
a surprise that the differential equation is the same, up to a map w — ► w, q — > q. The range 
of values for the problem of SM()12)I is | < q < < 0, while the sausage is characterized 
by ^ < cj < 1, iu > 0. The point is that Eq. 1)11)1 can be brought to Heun's form by the 
transformation £ = exp{2(y + u)}. The singularities are now located at {0,—1,—w, oo}, 
with w = exp{4it} and 24q = 6(1 + w) + k(u)(w — 1). The equation is actually the same 
as the one we find for the SM I)12|) . but the domain involved is the positive real line instead 
of the unit interval and the singularities are differently situated. By applying well-known 
transformation properties of Heun's equation (Ref. ^lEI!) we can reposition the domain 
on the unit interval (0, 1), the singular points being now 0,1, w = w/{w — l),oo and the 
new accessory parameter is given by 

(36) 



10—1 

Notice that the map {q,w} — ► {q,^} is involutory with w = oo as the only fixed point, the 
interval \ < q < 1 being mapped onto 1 < q < cxd. 

Hence we can use the same algorithm of the previous section in a different domain. The 
"magic" here is provided by the analyticity properties of the models involved. The complex 
shift transforming the "sausage" model into the SM 1)12)1 does not modify very much the 
eigenvalue equation, which turns out to be the same equation in a different domain. 

3. The Lame formulation 

The parameter A defined by Eq. I|3()j) is naturally related with a reformulation of Eq. 1)13)1 
close to Lame elliptic equation. If we define the modulus of Jacobi elliptic functions 

(37) k 2 = A = 1/(1 + exp(-4u)) 



then the substitution 

y _ u _ dn(z|fc 2 ) _ / sn(z|fc 2 )dn(zjfc2)" iTi 

{M) ~ ksn(z\k 2 ) ' ^ m "V cn(z|A; 2 ) ^ m 



/ a an \Z,Z\K ) TIL CO. \Z\Kr ) \ i 



maps the point y = oo to z = 0, the point y = — u to z = K, where K(k 2 ) is the real 
period of Jacobi elliptic functions, and it turns Eq. ((T3]) to the form: 

d 2 dn 2 (2z|£; 2 ) m 2 cn 2 (z|A; 2 ) 
dz 2 ~ sn 2 (2z|A; 2 ) + sn 2 {z\k 2 ) dn 2 (z\k 2 ) , 

with the boundary conditions ip m ~ z m 1/2 at z -> 0; Vm ~ (-K" - -z) 1/2 at z -> if. This 
equation can be studied analytically in two limits A: 2 — > (u — > — oo) and k 2 — > 1 (u — ► — oo). 
In the first case we can develop the standard perturbation theory near the exact solutions 
ipm n{z) = V sin 2z cos" 1 (z) Pn (cos 2z) where Pn a '^\x) are Jacobi polynomials. This 
perturbation theory gives the same IR expansion for the eigenvalues which was considered 
in previous section. 

In the opposite limit u — > oo, k — ► 1 and the real period K ~ — |log((l — fc 2 )/16) ~ 
2u+2 log 2 — ► oo. In this case the potential term in Ea. (|39|) is equal to m 2 almost everywhere 
and near the points z = 0; z = K it can be approximated with exponential in u accuracy 
by the potentials: 

(40) V{z) = \ — + m 2 coth 2 z 0<z-^K 

sinh 2z 

and 

(41) Vi(zi) = 2 h m2 tanh2 z i 0<zi = if-z<i<: 

sinh 2zi 

Both these potentials appeared in Ref . |27| . where the spectrum of CFT describing Witten's 
two-dimensional Euclidean black hole |2B] was studied. There, it was noted that potential 
V\(zi) is attractive and has the bound states solutions: 



(42) Vm,n = \/ tan hzi (coshzi) 2 ™ m+1 F (— n, -n + m, m - 2n; 1 - tanh 2 Zi) 

where F(a, b, c, z) is Gauss' hypergeometric function. These solutions are normalizable for 
integer n < (m — l)/2 and give the levels 

(43) \k = {m 2 - (2n + 1 - m) 2 \ n = 0, 1, [\m] - 1} . 

The corresponding eigenvalues of Eq. (j39|) approach these levels exponentially in u. 

The potential V{z) is repulsive and does not have normalizable solutions. For 2n > m — 1 
we parameterize K m „/6 = m 2 + p 2 n - Then solutions regular at z = and at Z\ = K — z = 
are found to be 

(44) 4>m,n{z) = (tanh z) m+ ^ (cosh z) ip F (|(1 + m - ip), \{l + m - ip),m+ l;tanh 2 z) 



(45) ip m ,n(zi) = (tanh £1)2 (coshzi) ip F (^(1 + m — ip),^(l — m — ip), 1; tanh 2 z\) 

Matching these solutions with the plane wave solution in the region <C z <C K we obtain 
the quantization condition: p„ = jir(2n — m+2)/(u+r m ) where r m = ijj(l) — ijj(^m+^), and 
ip(x) is the logarithmic derivative of the V function. This condition leads to the asymptotics 

(46) |sn(") =m 2 + 7T 2 {2U ~ ™ + §- + Q(l/n 5 ); n>±(m-l) 



(see Fig. |IJ). We note that for m^O this UV behavior is different from that for the sausage 
model eigenvalues, which is given by 

(47) \< 9 n(u) =m 2 + tt* ( " + 1)2 . 2 + 0(l/u 5 ); (n > 0) 

The sausage model eigenvalue problem of Eq. (jllj) can also be written in the elliptic form 
of Ref. Q31 with k 2 = 1 - exp(-4u): 

( d 2 cn 2 (2#g) , m 2 dn 2 (z|fc 2 ) \ , ^ n k 2 s , 

1 j V dz 2 sn 2 (2z|fc 2 ) + sn 2 (z|/t 2 )cn 2 (z|fc 2 )y'^ ?n ' n " 6 ^ m '" 

with boundary conditions -0m ~ z m+1 / 2 at z — ► 0; i\) m ~ (if — ^) m+1 / 2 at 2 — ► if . In the 
UV limit u — > oo the potential term in this equation at the both ends tends to V given by 
Eq.pnj). The normalizable solutions do not appear and both asymptotics can be described 
by the eigenfunction (|44|). The quantization condition for parameter p n in this case leads 
to the asymptotics Eq.(|47|). 

4. Numerical analysis 

The matrix representation introduced in a previous section, while useful from the analytic 
viewpoint, is not the best choice if we want to compute the spectrum beyond perturbation 
theory. Actually an n-dimensional truncation of the matrix given in Eq. 1)22(1 is going 
to be essentially equivalent to n-th order perturbation theory. We are now introducing 
another transformation of Eq. ()13|) which will allow us to efficiently explore the whole range 
— oo < u < oo. 

We start from the fact that in the limit u — * — oo (and m = 0) there is a simple solution 
with 



(49) = VT - e- 2 ^+«) = a(y) , 

with k(u) — > 6. Hence it seems promising to look for a solution of the form 

(50) vP(y) = <x(y)V>(2/) 

To find the new differential equation it is convenient to write down the functional (7i) whose 
critical points are the eigenvalues: by denoting x(u) = 1 — sinh2u sinh2y and recalling 
p(y) from Eq. we have 

f°° (^(m 2 - d 2 )^ - XP~ 2 ) ^ 2 ) dy 

( 51 ) W = -1*2 / 

L u p * d y 

Now, by inserting Eq. I|50|). after an integration by part, we find 

_ I-u [* 2 r + K - X<? 2 /P 2 - (™y + m 2 a 2 ) ^] dy 

The variational equation resulting from this functional 

(53) — 7— (a 2 -r- 1 + = cosh 2-u <r 2 ib 

dy \ dy J 

suggests to introduce a new coordinate £ such that dy = u 2 d£; we find 

(54) y = -u + ±log(l + e 2 «) 
and the equation simplifies to 

(55) + a 2 V(£, u)^ = \k{u) W(£, u) ^) 
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with 

m * - e2{((i+e4, ];:^,;T^:r +3e2 "' + ^ 2(i+t - h « 2 

e 2 Hl + e 4u ) 



(5?) W (l + e 2 «)(l + e 4 « + e 2 €) 
The crucial property of this new formulation is that it is now regular on the whole real axis. 
With appropriate boundary conditions (Neumann's b.c.)^ the eigenvalue equation can now 
be easily solved by standard sparse-matrix generalized eigenvalue routines ** 





Figure 2. The generalized potentials V, W (upper) and the hrst eigenfunc- 
tions ifi 2 for u = 10, m = 0. 

See Fig.s Q2I3JI for some typical waveforms. The density includes the measure appropriate 
for the new variable £, namely er(£) /p(£)- In the limit u — ► — oo we easily recover the 
discrete spectrum 

(58) \k = {(2j + l) 2 - m 2 | j = n + ±m, n = 0, 1, 2, . . .} 

as we already know from perturbation theory (see Eq. (jSHJ)- in the other limit, u ^> 0, 
the spectrum can be described as a monotonously decreasing flow toward k/Q = m 2 , except 
for a finite number of eigenvalues which have a value less than m 2 : these peculiar "bound 
states" are given by Eq. (|43jl and they are easily reproduced numerically (see Fig. 0J. We 
easily check that they nicely agree with the asymptotic formulas already given in Eq. (|47|). 



5. Matching to tba data 
As it was already mentioned in Introduction for the special values of parameter 

^ £=^ ; Ar = 3 - 4 -" 

"it is somewhat tedious to trace the b.c. back from the original equation; suffices to say that, intuitively 
speaking, since the constant solution is exact in the limit u —> —oo, Neumann b.c. are the natural ones. 
**We successfully used the routine eigs in Matlab. 
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a) 




-10 -5 5 10 15 20 



Figure 3. The densities at u= 10, m = 10: a) "bound states", b) the "continuum". 



to 
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II 
LU 




-5 5 10 15 20 

U 



Figure 4. The asymptotic behavior in the UV for the low lying states at 
m = 0: lines — > Eq. i)46'|) . circles from numerical solution Eq. i)55|) . 



the effective central charge of the SM (|12j) admits the exact (i.e. to all loops) calculation 
using the tba equations. These equations can be derived from the factorized scattering 
theory of right (r) and left (I) moving massless particles which form the spectrum of our 
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U 



Figure 5. The spectrum How for m = 10; dashed line shows the continuum threshold. 



MR cosh [5 




1 N-l 

Figure 6. The extended Dynkin diagram for the Sausage model. 



SM. The factorized scattering theory of massless particles is characterized by two-particle 
S-matrices S rr (f3), Su(f3), and S r i{j3) where (3 is the relative rapidity of scattering particles. 
We suppose to discuss in more details this scattering theory (as well as perturbed CFT 
approach to our SM) in another publication. Here we only note that each of these three S- 
matrices coincides formally with the S'-matrix for the massive particles in Zfq parafermionic 
CFT perturbed by the parafermionic operators. This S'-matrix is described in details in 
Refs.j22l I30j . This scattering theory results in TBA equations which form the common 
system of N + 1 coupled non-linear integral equations for N + 1 functions e a (/?) of rapidity 
variable — oo < j3 < oo. The tba system has the form: 

N 

(60) pa ( (3 )=e a + — / ^^ a6 (/3-/3 / )log(l + e- £ ^))d/3 / ) 

n J 6=0 
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1 N-l 
Figure 7. The extended Dynkin diagram for the SM lfEZl) 



where (f a b(P) = cosh/3, T being the incidence matrix of the extended affine Dn 

Dynkin diagram (see Fig. 7) and the source term 

(61) p a (P) = ±RMeMP)SaO + ^RMeM-P)$aN. 
The effective central charge can be calculated as 

(62) c(R) = \ f Paid) log(l + exp(-e a (/3)) d/3. 

The incidence matrix of these tba equations is similar to that for the sausage model (see 
Fig. 6), the only difference coming from the source terms. For the massive sausage model 
(without topological term) p a {f3) = RM cosh(f3) 5 a o- 

The UV behavior of c(R) is determined only by the structure of the incidence matrix 
and we can find that in both cases effective central charge approaches to the limiting value 
c uv = c(0) = 2 logarithmically (see Ref.^^) in agreement with Eqs. (|46|47[) . The analysis 
of TBA equations with source p a given by Eq, (|61|) shows that in the IR limit cm = c(oo) = 
2 — 6/ (N + 2) coincides with central charge of parafermionic CFT and the IR corrections to 
this value have a structure: 

e { R,N) _ 2- — O 2 (A0(^_) + 

N + 2 \ 12/(JV+2) 



(63) ^{-mr) 

The IR asymptotics of our theory can be described by the methods of perturbed CFT. This 
field theory is characterized by integrable perturbative operator which belongs to the space 
of fields of the parafermionic CFT and has the dimension A per t = 1 + 2/(N + 2). It gives 
us the possibility to calculate analytically the first corrections to the expansion ()63|) . The 
exact values for the coefficients b2(N) and 63 (iV) are presented in Appendix C. 

The effective central charge c(R, N) was computed using the tba equations numerically 
for several values of N (N = 5, 7, 11, 15, 23, ...). At large N the central charge is predicted 
to be given by 

(64) c(R, N) = 2 - + 0(1/N log N) 

To verify this we have to relate u to the parameter MR entering the tba equations. Within 
one loop approximation we have a freedom which can be used to fit the data in the best 
way. Asymptotically we expect u ~ \og(N/MR)/N, but finite N corrections are present 
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and may be important. Empirically we find that a rather accurate choice is the following 

(65) u = N^ 1 log{N/MR) 

(66) N cS = \J (N + 2)(N - 2 tanh(4 log(N/MR))) 

which is used to build the plot of Fig. ED It is quite evident that the data are increasingly 
well matched by k as N increases (in the deep UV or IR the agreement is even better). 




nft«nV9f> 6<<^& ■ iga O BjD (0 nDn CD BP O 
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u=log(N/MR)/N 
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Figure 8. Matching k to the tba data. 
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Figure 9. A close-up view of Fig.9. 
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Finally let us note that the exact values of the coefficients 62 (N) and 63 (N), given in 
Appendix C, have been reproduced with a high degree of accuracy by fitting the TBA data 
at various values of N. 

Conclusions 

We have found a unified treatment of the general spectral problem for the sausage model 
and its variant SM ()12() . The two formulations correspond to the same Heun equation 
defined on different domains in complex plane and they are linked by a projective trans- 
formation. The spectral function k(u) can be studied both numerically by diagonalizing 
a discretized form of the differential operator, or analytically by a power series expansion 
around the IR point. The perturbative expansion can be pushed to high orders and it turns 
out to be convergent in the whole physical domain. Its asymptotic behavior at high order 
is compatible with the leading UV behavior, which has been computed analytically. In the 
ultraviolet regime, the spectrum, besides the expected continuous component, contains a 
set of bound states with angular momentum higher than one. The scaling functions for the 
ground state has been compared to the central charge computed via tba equations: a very 
good agreement was found, adding good evidence for the interpretation of the SM (|12j) (at 
v = 4tt/N) as the field theory describing the RG flow to the Zn parafermionic CFT. 
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Appendix A 

We recall some standard results from Raileigh-Schroedinger perturbation theory for non- 
degenerate levels, which should be found in any textbook on Quantum Mechanics. Sta- 
tionary perturbation theory is most efficiently formulated as an iterative algorithm. Let 
H = Hq + e V, Eq any one of the unperturbed eigenvalues, with eigenstate \Eq). Then the 
perturbed eigenvalue E(e) can be expanded as 

(67) E(s)=E + J2tn£ n 

n>l 

while the eigenstate is given by a vector series in terms of auxiliary vectors \rj n ) as 

(68) \E(e)) = \E ) + J2e n \r ]n ) . 

n>l 

The expansion coefficients can be computed through the following recursive algorithm: 



Let 1 770) = l-Eb) and TZq = — —J — °^ ; then for any n > 1 we have 

Hq — En 



(69) 




{E \V\rin-l) , 

= TZq {— V|77 n _i) + Y2=l $k\rin-k)} 



The algorithm can be easily translated in any symbolic manipulation language. In prin- 
ciple we should take care of defining the matrix V in an orthonormal basis; this would 
introduce some square roots in our matrix Eq. (|22|) . while at the end the coefficients turn 
out to be rational. Actually the following lemma tells us that we may comfortably work 
with the unnormalized basis. 

Lemma: Let H, Hq, V,\Eq), E(e),TZ, 5 n ,\ri n ) be defined as above. Let us assume further, to 
avoid any convergence problem, that V be a finite-band matrix along the main diagonal. 
Let S be any nonsingular Hermitian operator commuting with Hq. Then we may substitute 
V with SV S 1-1 in the iterative algorithm Q69J) . leaving everything else unchanged. 
Proof: Let 5' n , \r]' n ) denote the sequence constructed by inserting SVS -1 into eq. (|6^|) . Since 
S commutes with Hq, we must have S\Eq) = sq\Eq) for some non-zero sq. Hence 

(70) 5' n = {EQ\VsQS- l \r,' n _ l ) . 

Let us define |r^) = so5' _1 | r /n-i)- eas y *° check that the sequence {5' n ,\r]")} satisfies 
the same recursion as {5 n ,\r]n)}, and moreover |t/q) = sqS~ 1 \Eq) = \Eq), hence the two 
sequences must coincide. 

Note: the statement in the lemma is strictly perturbative. H and S H S^ 1 could be in- 
equivalent as operators, since S and/or S^ 1 may be unbounded, still they share the same 
perturbative expansion. 



Appendix B 

We report the explicit expression of the series expansion of K mj n for small values of (n, m). 
Of course this is just a sample; the code can generate them to any order, the only limitation 
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^0,0 
«1,0 = 



being computer's physical memory and time. 

g \2 1 \3 229 \4 109 \5 62999 \ 6 _ 20159 \ 7 _ 299803787 \ 8 
u A 2 A 720 A 480 A 362880 A 145152 A 2612736000 
72503387 \9 _ 173336436487 \10 , nf^ 1 ") 
746496000 A 2069286912000 A t <^ A ) 
18_4\_16\2_352 \3_ 1972 \4 _ 17408 \5 _ 701314 \6 
10 9 A 405 A 3645 A 45927 A 2460375 A 

34835788 \7 _ 204567413 x 8 _ 1588447666493 x9 
155003625 A 1116026100 A 10358117240625 A 
4782354354298021 \10 , nAll\ 
36543437624925000 A / 
— — 94 X _ 27 \2 _ 27 \3 _ 23949 \4 _ 34047 \5 _ 370287 \6 

ft o, l — 01 A 5 10 A 14000 A 28000 A 400000 A 

826209 \7 _ 146655243891 \8 _ 35351959491 \9 
1120000 A 241472000000 A 68992000000 A 
197594782006203 \10 , nfX 11 ) 
448448000000000 " r \ A > 
qn — X — 9664 x2 _ 7627904 \3 _ 217386688 \4 _ 173655964928 \5 
au 5 A 1125 A 1771875 A 79734375 89701171875 A 

485256409132928 \6 _ 955858372577612032 \7 _ 176847414696606187696 \8 
329651806640625 A 815888221435546875 A 183574849822998046875 A 

12815580494902423265411456 \9 i /O('X 10 'l 
15786519210528717041015625 A / 



Appendix C 

In this appendix we give the exact values for the first IR correction to the levels e m j(R) 
and two first corrections to the effective central charge or to the ground state energy eo(R). 
We express these corrections in terms of parameter M entering the tba equations. These 
corrections can be calculated analytically using the methods of integrable perturbed CFT. 
We consider the case v/4tt = 1/N. In this case the IR limit of SM (|12|) is described 
by parafermionic CFT. The conformal dimensions of the primary fields in this CFT are 
characterized by two quantum numbers m and j (j = \m/2\,\m/2\ +1... < N/2) and have 
a form: 

(71) A mj = 3{J + 1} - ^ 

V ' m) JV + 2 W 

It is convenient to introduce D m j = (N + 2)A m j. As it was noticed in the Introduction, 
the limiting values of (e m j — eo)/24 should coincide with A m j. With the first IR correction 
these values are: 

(N + 2)(e mi (R)-e (R)) n2 h(j,N) (N + 2\^ N+2 ^ 

{U) 6 ~ mj mj j(j + l)\MRj + - 

where the coefficient bi(j, N) can be expressed through the function g(x) = T{l+x) /Y{l—x) 
and has a form 

(73) h(j,N) = ^ 2 g (^) 2 g (|g) 4/(^2)^ 

(^ + 2 ) 2 9(7^2) g(wk) 

We note that in the large limit the left hand side of Eg. 1)72(1 tends to (K m j — k,q)/6. For 
4:D m j and b\ (j, N) we have 

(74) 4D mj = (4j(j + l)-m 2 )(l + 0(l/N)), h(j, N) = 2 + 0(1/N) ; 

if we now define the parameter u by the relation: ( j^ 2 ) 4 ^^ 2 ^ = exp(4u) = A/(l — A) we 
can see that the first term in the expansion given in Eq. (|33|) coincides with Ea.(|72jl at one 
loop accuracy. 
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For the effective central charge c(R, N) (or for ground state level eo(-R)) it is possible to 
calculate analytically two further IR corrections. Namely, the coefficients 62 and 63 in the 
expansion are 

N 2 (N - 2) 2 0(1) g(Jrn) (Svr) 8 /^ 2 ) 

<™> * + - J1 Z:Z f 1 T' = -f + 

2 (JV + 4) 4 (JV + 8) 2 5(77+2) 5(7v+2) 2 
It is easy to see from these equations that with one loop accuracy the IR expansion for the 
function kq(u) = 6 — A 2 — A 3 /2 + ... coincides with exact IR expansion for the function 
(N + 2)(2-c(R,N)). 

In the UV limit the leading term for the levels e m j(R) (for j > m — 1/2) can also be 
calculated exactly and has a form: 

(77) Ne mj (R) = ^ + ' b" - " + ! W - 2 > + 0(1/3 



2Z£(fl) 



where Z m (R) = log(8vr(7V - 2)/RM) + (N - 2)(^(1) - VK 1 ^)) + V'(l)- This asymptotic 
behavior coincides at one loop accuracy with Ea.(|46|). where n = j — m/2. 
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